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We propose an integrable system of coupled nonlinear Schrodinger equations with cubic-quintic 
terms describing the effects of quintic nonlinearity on the ultra-short optical soliton pulse propagation 
in non-Kerr media. Lax pair, conserved quantities and exact soliton solutions for the proposed 
integrable model are given. Explicit form of two-solitons are used to study soliton interaction showing 
many intriguing features including inelastic (shape changing) scattering. Another novel system of 
coupled equations with fifth-degree nonlinearity is derived, which represents vector generalization 
of the known chiral-soliton bearing system. 



I. INTRODUCTION 



Optical solitons have promising potential to become principal carriers in telecommunication due to their capability 
I ■ of propagating long distances without attenuation and changing their shapes ■ Therefore considerable attention 
is being paid theoretically and experimentally to analyse the dynamics of optical solitons in optical waveguides 
(for example silica fibers) under different contexts Qd^]- Such investigations are helpful for realizing optical soliton 
, applications, particularly in soliton-based optical communication systems j^] and nonlinear optical switches The 
waveguides used in such optical systems are usually of Kerr type Q . Consequently the dynamics of light pulses are 
described by the nonlinear Schrodinger (NLS) family of equations with cubic nonlinear terms However, as the 

intensity of the incident light field becomes stronger, non-Kerr nonlinearity effect comes into play and due to this 
| additional effect, the physical features and the stability of NLS soliton can change ||. 

• The way through which non-Kerr nonlinearity influences NLS soliton propagation is described by the NLS family 
of equations with higher-degree nonlinear terms |p|-|l6| . Therefore investigations on these evolution equations become 
important from a theoretical point of view. Particularly this importance has received a boost after the experimental 
* r j H ', observation of multistability of solitons in non-Kerr fiber |l7]] . In general the models proposed in the literature [PU9j-|T^1 
for describing the non-Kerr effects are not completely integrable and can not be solved exactly by the inverse scattering 
transform method. In such non-integrable systems, therefore, the details of soliton-interaction during collision can 
not be described exactly and hence are still open to debate. However, numerical stimulations |Q show that even the 
slightest change from the Kerr nonlinearity results in the two solitons annihilating each other, merging or creating 
many new solitons, depending on the initial inclination of the two solitons and their shapes. But besides the important 
problem of computer time, the numerical approach is not very appealing in the sense that it is not a simple task to 
get physical insight from purely numerical experiments. The idea therefore is to use approximate analytical methods 
such as perturbation technique, variational method, etc., in order to compensate for the lack of exact results p9| . By 
treating the quintic nonlinear terms due to non-Kerr effect as perturbations of the cubic NLS equations, i.e. restricting 
the effects of quintic nonlinearity to be less predominant than the cubic terms, the NLS equations are studied both 
analytically and numerically in j|,|nj|l3],[l6| ■ 

In this paper we have obtained for the first time (to our knowledge) an integrable system of coupled NLS equations 
including cubic-quintic terms describing the effects of quintic nonlinearity with arbitrary coupling, which generalizes 
the coupled hybrid NLS equations with cubic nonlinearity |2(],^l| . Lax pair as well as infinite set of conserved quantities 
are derived for the proposed integrable model. We also find exact soliton solutions for our model and using the explicit 
form of the two-soliton solution we study the associated soliton collision. A remarkable interrelation between the newly 



proposed integrable model and the celebrated Manakov model 22 helps to use the recently discovered general two- 
soliton solution of the later model (2^] . This reveals the fascinating occurrence of shape changing inelastic soliton 
collision also in the present model, in addition to some other interesting features. We believe that such a study 
using higher-order solitons (multisolitons) becomes important in the light of the proposal by Hasegawa and Nyu |2J] 
regarding " eigenvalue communication" , in which the information may be transmitted by higher-order solitons and a 
recent one sugested by Jakubowski, Steiglitz and Squiev ]2^| on nontrivial information transmission system, which 
uses the more general two-soliton solution of the Manakov system j22| derived by Radhakrishnan, Lakshmanan and 
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Hietarinta in |23|] . We also find that the Hamiltonian of the present integrable model is associated with a non-canonical 
Poisson bracket (PB) structure. However, using the same Hamiltonian with canonical PB relations we derive another 
coupled system with quintic nonlinearirity, which may be transformed to a vector generalization of the chiral-solitonic 
model of Aglietti et al |53| . 

The plan of our paper is as follows. The basic evolution equation with cubic-quintic nonlinearity which describes 
the soliton evolution in non-Kerr media with parabolic nonlinearity is discussed in Sec- II. Sections III is devoted to the 
proposed model and presents its integrability property by explicit construction of the Lax pair as well as the hierarchy 
of conserved quantities generated through a recurrence relation, which in turn is derived from a coupled Ricatti 
equation using the Lax operator. Sec IV gives the exact soliton solution for our cubic-quintic nonlinear evolution 
equation and using the explicit analytic solution we study here the collision process of solitons so as to understand 
the influence of quintic nonlinear effects on the Manakov model. Section V establishes the interrelation and gauge 
transformation between our system and the Manakov model including the anyon-like nonultralocal PB for our model. 
Sec VI presents a vector generalization of the chiral soliton bearing system starting from the same Hamiltonian and 
discusses some of its interesting features. Sec VII gives concluding remarks. 



II. BASIC EVOLUTION EQUATIONS 



Generally when high optical intensities (or materials with high nonlinear coefficients even at moderate optical inten- 
sities, for example semiconductor doped glasses, organic polymers, thin liquid filled capillaries, etc.) are considered, 
it is necessary to take into account higher power nonlinearities arising from an expansion of the refractive index in 
powers of intensity / of the light pulse : n = no -\- n^I -\- n^I 2 + where no is the linear refractive index coefficient 
and 712,^,4..., are nonlinear refractive index coefficients [p|]ic|]. 

In the case of n = Uq + n%I + n^I 2 , the wave equation for high-intensity light pulse propagation in an isotropic 
single-mode optical fiber with a circular cross-section and fiber axis z can be written as 
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where c is the speed of light, the linear part Dl and the nonlinear part Dnl of the electric field displacements are 
related to the electric field E(r,t) by the relation Dl = J °° e(t')E(t — t')dt' and Dnl = ^2\E\ 2 + e^E^, in which 
e = n 2 , £2 = 2n2no and e4 = 2rnriQ. 

A solution of Eq.(I) is sought in the form 



E = eR(r)A{z,t)e 



i(3z—iujt 



(2) 



where e is a unit vector in the direction of wave polarisation, R(r) describes the transverse field modes, in which r is 
a two dimensional vector in the x-y plane and A(z, t) is a slowly varying amplitude. Here we assume that -R(r), which 
is mainly defined by the linear effects, corresponds to the modal distribution of the fundamental fiber mode HEn, 
for simplicity. Then from the equations (I) and (2), assuming the temporal dispersion of the dielectric permitivity 
to be small, using the slowly varying envelope approximation and following the procedure in section 2.1 of the 
follwing nonlinear partial differential equation for A(z, t) can be obtained, 
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where the subscript to of the wave number k ( i.e. k u ) indicates differentiation of k with respect to to, the subscripts t 
and z of A indicate differentiation of A with respect to the coordinates t and z respectively and the numerical values 
of the parameters ao and /?o depend on the form of the function R(r). 

It is convenient now to transform the above equation to a reference frame moving with group velocity v g , and 
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nonlinear propeties of the fiber and \Aq\ is a measure of the maximum amplitude of the input pulse. Now equation 
(3) takes the form 



iq z + qtt + 2\q\ q + j\q\ q + ijiq m + «7a(M ?)t + *73(M Q)t = 0. (4) 

Equation (4) describes the effects of quintic nonlinear terms proportional to the real parameters 7 and 73 on the 
dynamics of pulse envelope allowing self-phase modulation and higher-order linear and nonlinear dispersions. For 
pulse widths greater than 100 fs, one can neglect the last three terms of (4) and the resulting equation is a well- 
studied simple normalized NLS equation with cubic-quintic nonlinear terms. Further the system (4) is a special 
case of the dynamical equation considered for the fiber system with saturating nonlinearity p2| |. 

It is of further interest to extend the above analysis to include multimode effects. For this purpose, there are several 
ways to generalize equation (4) to a set of coupled equations depending on the physical situations. A fairly general 
form of coupled nonlinear Schrodinger (CNLS) equation with cubic-quintic nonlinearity is 

iqu + qitt + 2{\qi\ 2 + B\q 2 \ 2 ) qi + 7(ki| 2 + B\q 2 \ 2 ) 2 qi + pq x + nq 2 
-i^htqutt +72[(|<?i| 2 + B\q 2 \ 2 )q lt + (qiq\ t + Bq 2 q\ t )q x 
WiQuQi + Bq 2t q* 2 )qi}+l 3 [(\qi\ 2 + B\q 2 | 2 )V] t ] = 0. 
iq2z + q2tt + 2(B|<7i| 2 + \q 2 \ 2 )q 2 + 7(-B|<?i| 2 + \q 2 \ 2 ) 2 q 2 - pq 2 + nqi 
-ilAliqmt +7 2 [(B|<7i| 2 + \q2\ 2 )q2t + {Bqiq* lt + q 2 q 2t )q 2 

+i{Bq uq l + q2t q*)q 2 ]+ l3 [{B\ qi \ 2 + \q 2 \ 2 ) 2 q 2 ]t] = 0. (5) 

A nonlinear direction coupler with quintic nonlinearity (or parabolic nonlinearity coupler) has B = p = p, = [ pc[ . 
For p, = 7 = 0, Eq.(5) acts as a mathematical model for a periodically twisted elliptical birefringent fiber ]27|] . If 
7 = p = K = 7i=73 = 0, and 7' =1 then Eq.(3) becomes the coupled hybrid nonlinear Schrodinger equation |2(]] 
used to investigate the effects of birefringence on pulse propagation in the femto-second range. In the absence of 
quintic nonlinear terms proportional to the real parameters 7 and 73, soliton interaction supported by system (5) has 
been studied by deriving higher-order soliton solutions under the parameteric restrictions B = 1 and 371 = 72 [ p8[ . 
One may also note that when B = 1 and 7' = 1 the linear coupling terms proportional to the parameter p and k in 
Eq.(5) can be removed without affecting the other terms by using the transformation 

qi -> cos(^)e irz q 1 - sin(^)e- irz q 2 , 
6 9 

q 2 -► sin(-)e irz qi + cos(-)e- irz q 2 , (6) 

where T = (p 2 + n 2 )i and 6 = tan^ 1 ^). If the nonlinearity is restricted only to cubic terms corresponding to pulse 
widths greater than lOOfs, one obtains the celebrated integrable Manakov model p2| , |35| 

iqiMz + q\Mtt + 2(\qiM I 2 + |<?2m| 2 )<Zim = 0, 

iq2Mz + q2Mtt + 2(|<jia-/ | 2 + |<72M| 2 )<?2M = 0. (7) 

There is a large amount of theoretical work devoted to the CNLS family of equations with cubic nonlinearity. 
However, to the knowledge of the authors, to date the CNLS equations with non-Kerr nonlinearity have received very 
little attention in the literature, particularly in connection with the integrability aspects. In the following sections, 
by identifying one such integrable nonlinear evolution equation, we derive the two-soliton solution so as to get some 
idea about soliton interaction in non-Kerr media. 



III. INTEGRABILITY PROPERTY OF THE PROPOSED MODEL: LAX PAIR AND CONSERVED 

QUANTITIES 

It is evident that eq. (5) does not exhibit the explicit rotational symmetry in the internal space spanned by the 
vector (gi,g 2 ). However, for B = 1 such a symmetry is restored. Assuming further that 71 = 73 = 7' = 0, (5) can be 
reduced to the following quintic generalization of the coupled cubic NLS equation, 
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iqiz + qut + 2(|?i| 2 + 1^2 | 2 )^i +7(l?i| 2 + \i2\ 2 ) 2 qi + pqi + Kq 2 - ^72 

[iki\ 2 + l92| 2 )<Zit + (qiqtt + Q2Q2t)Qi] = °> 
iq2z + qitt + 2(\qi\ 2 + \q2\ 2 )qi + 7(l<7i| 2 + lQ2| 2 ) 2 g2 - pq 2 + Kq\ - i/172 

[dftl 2 + k 2 | 2 )g 2t + {q\q* t +g2<Z2t)<Z2] = 0. (8) 

The terms p and k can be removed from (g) by using transformation (6). Eq. (^) without quintic nonlinearity was 
investigated in pi]. However the remarkable fact is that (|h itself can be shown to be exactly integrable. Our proposed 
model is a further generalization of (||) and naturally of (0), where the internal rotational symmetry is broken again 
and more parameters are introduced with arbitrary values, which can be chosen conveniently to suit the real situations. 
The model can be given as 

iqiz + qut + 2(\qi\ 2 + \q 2 f)qi + (piki| 2 + P2 1^2 | 2 ) 2 Qi + 2p 2 [(n - Pi)|gi| 2 + (t 2 - p 2 ) 
I92I 2 ] \q2\ 2 qi - 2i[(p 1 \q 1 \ 2 + p 2 \q 2 \ 2 )q 1 ] t + %{pxq\qu + p 2 q 2 q2t)qi = 0, 
iq2z + qm + 2{\qi\ 2 + |?2| 2 )?2 + (n|<?i| 2 + T 2 |g 2 | 2 ) 2 <7 2 + %n[(pi - ti ) l^i | 2 + (p 2 - r 2 ) 

I92I 2 ] ki| 2 <72 - 2i[(r 1 \q 1 \ 2 + T 2 \q 2 \ 2 )q 2 ] t + 2i(r 1 qlq lt + r 2 q 2 q 2t )q 2 = 0. (9) 

where pi,p 2 ,Ti and t 2 are real free parameters. It is evident that with a symmetric reduction p\ = p 2 = n = r 2 , 
we can recover (0) from (^|), while a different reduction with qi ~ q and q 2 = (or gi = and q 2 = q) yields the 
integrable Kundu-Eckhaus equation 

iq z + qtt + 2\q\ 2 q + p 2 \q\ 4 q - 2i Pl (\q\ 2 ) t q = 0. (10) 

Importantly this generalized model @ turns out also to be exactly integrable. For establishing the integrability 
property of the proposed system, which consequently proves also the integrability of the reduced model (^), we find 
the Lax pair (L, M) associated with (||) as 

(-iX qi q 2 

-q\ -i9 u + iX | , (11a) 
~q* 2 -iflat + *A 
/ [-2iX 2 + i{\ qi \ 2 + \q 2 \ 2 )], 2X qi + iq lt + 6 ltqil 2Xq 2 + iq 2t + 9 2t q 2 
M=\ -2Xql+iq* lt -6 lt ql [2iX 2 - i\qi\ 2 — i9\ z ], ~iq\q 2 ). (lib) 

\ -2Xq\ +iq^ t -9 2t q2, -%?2> [2iX 2 - i\q 2 \ 2 - i9 2z ] 

where 

01= / (pi|<Zl| 2 +P2|<72| 2 X 9 2 = [ (nlq^ 2 + T 2 \q 2 \ 2 )dt' . (11c) 

J — oo J — oc 

Here A is the spectral parameter. It may be easily checked that the zero curvature condition L z — M t + [L, M] = 0, 
with the explicit Lax operators (11), yields Eq. ([)]). In section V we give another evidence of its integrability by 
relating system (||) to the integrable Manakov model through a gauge transformation of the pair (11) to the Manakov 
Lax operators. 

Integrable systems, as is well known, possess infinite number of conserved quantities in involution, of which usually 
the lower ones are of physical importance. Explicit forms of such conserved quantities for the integrable system (j^) 
can be derived from a recurrence relation obtained from the Ricatti equation related to the Lax operator. For this 
purpose we use the linear system related to (11a) 

* t (A,t)=L(A,t)*(A,i), $ = (0i,02,0 3 ) (12) 

and observe that 

In o(A) =ln0 ie - lA Voo =J2 C « A ~" 

n 

serves as the generator of the conserved quantities {c„} through an expansion in the spectral parameter A. The first 
equation of the system (O) thus yields the relation 
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dt( qi rp + g 2 rW), n>l (13) 

-co 

with the expansion = J2n=i r£ o) A _n , a = 1, 2, where we have denoted r (1) = and r (2) = |a. For finding now 
the infinite set of conserved quantities we may use the rest of the equations of (12) to derive a set of two coupled 
Ricatti equations for rW(A) and r^ 2 ^(A). Expanding in powers of A as mentioned above we obtain the recurrent 
relations 

n—1 n—1 



MT& - r« + ^rw + qi j: r&if + , 2 £ r« rf >, (i4») 

n — 1 n — 1 

-ar^ = r£ + i0 2t if ) + q , £ iW + 9i E W ■ ( 14b ) 



i=l i=l 



with rj '' = — This gives finally the conserved quantities in the explicit form as 



1 



Cl = --/ dt(\ qi y + \q 2 n (15) 



X 



C2 = -J / dt[-l(q 1 q* lt +q 2 q; t )+p 1 \q 1 \ 4 + T 2 \q 2 \ A + (p 2 +T 1 )\q 1 \ 2 \q 2 \ 2 } (16) 

i f +ao 

c 3 = - g / + «M«t) + (kil 2 + M 2 ) 2 + 

l {\qi\ 2 N lt + \q 2 \ 2 N 2t ) + 2l(N iqi q* lt + N 2q2 q* 2t ) - (N 2 \ qi \ 2 + N 2 \q 2 \ 2 )}, (17) 



etc., where Ni = 6 U = pi\ q \\ 2 + p 2 \ q2 \ 2 and N 2 = 9 2t = t"i | (?i | 2 + t- 2 | (72 1 ■ The above conserved quantities in ( |l5f|l7| ) may 
be interpreted in terms of the number operator N, the total momentum P and the total energy or the Hamiltonian 
of the system H. However it is intriguing to remark that since here the fields q a ,qa do not have canonical Poisson 
bracket relations some care has to be taken in deriving the equation of motion (j^) from the Hamiltonian (j^) . In 
fact the Poisson bracket structures of the fields, which are derived in sec V, show an interesting anyon-like feature. 



IV. EXACT SOLITON SOLUTIONS AND SCATTERING OF SOLITONS IN THE GENERALIZED 

MANAKOV MODEL IN NON-KERR MEDIA 

The proposed system (0) with quintic nonlinearity allows exact ./V-soliton solutions, which can be found for example 
by Hirota's method following the same procedure as in the Manakov model. However a more direct and convenient 
way is to use the known solutions of the Manakov model (0) themselves for constructing the soliton solutions of (||). 
This is possible due to the interrelation between these two models, which will be established in the next section. Thus 
we find the explicit 1-soliton solution of (^) in the form 

( quq2 ) = (ae i5ltanh ( l/ ('- , ' 2+ ' 5 »,^e i52tanh(,y(t -^ +5)) )sec/i(i/(i- vz + 8))e i{Kt+u>z) , (18) 

where different parameters of the solution arc related with the spectral parameter A = v + in and the parameters of 
the model as 

v = 2k, uj = v 2 -k 2 , 8 1 = l( Pl \ a \ 2 + p 2 \{3\ 2 ), S 2 = i(n \a\ 2 + r 2 \{3\ 2 ) 

together with a constant phase 8. Comparing with the Manakov soliton we see that there is an interesting phase 
change in the carrier wave. The plane wave like character in the Manakov model has been deformed into a wave 
suffering compression and rarification of the phases in a kink-like profile (see Fig. la,b). This shows that the effect of 
quintic nonlinearity of our model appears in the soliton phases and therefore in the derivatives of the soliton profile 
9aa;j 9at i which must also change the momentum and energy of the soliton. 

Now exploiting the higher-soliton solutions of the Manakov model it is also possible to find higher-solitons for (jj^) 
in explicit form. The Manakov model (fjj) has received considerable attention in recent years |^,^,^3| in order to 
understand the soliton collision in birerefringent fiber. However the importance of finding higher-soliton solutions in 
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explicit form has been understood only quite recently 1 23 |25l] . By constructing the most general two-soliton solution 



of the integrable Manakov model, two of the authors (R.R and M.L) and Hietarinta j23| have shown that the soliton 
in birefringent fiber can in general change its shape after interaction due to an intensity redistribution among the 
modes, eventhough the total intensity remains conserved. This shape changing collision arising essentially due to the 
change in polarisation angle helps to realise the exciting possibility of switching between components. (However, the 
standard shape-preserving collision property of the (1+1) dimensional soliton system is recoverd, when restrictions are 
imposed on some of the free parameters in the two-soliton solution |^3| ) . Recently using the shape changing collision 
concept, Jakubowski, Steiglitz and Squier |2^] have designed sequences of solitons operating on other sequences of 
solitons that effect logic operations and suggested nontrivial information transformation system. 

Now in order to investigate the implication of this property of soliton when the additional cubic and quintic nonlinear 
terms are included, we construct the two-soliton solution of the system (^|) and expect again a non-trivial change in 
the soliton phase as in the case of 1-soliton solution. Such a two-soliton solution expressed compactly through that 
of the Manakov system (Q) (qiM,Q2M) takes the form (see eq.(24) below) 

„ _ „ J f (pi\qiM\ 2 +P2\q2M\ 2 )dt _ i f (ti | go,/ 1 2 +r 2 1 92M | 2 )d* / 1Q \ 

qi — qiMe ■> , q2 — q-iuz J > (tyj 

where a general two-soliton solution p3| of the Manakov model is given by 

{qiM,q2M) — ^ _|_ em+Vl , +Rl + em+rl2 , +So + eVl * +rl2+So * _|_ e r, 2 + V2 *+R 2 _|_ e i)i+t)l*+ij2+**+fl3 ' ^ ' 

in which m = kj(t + ikjz) , j = 1,2, e s ° = = ^ , = = (kl+ fcfe +ha) (aiK 21 - a a « u ), 

e S2 = {k2+ %H$\ +kV (<*2Ki 2 -o llK22 ), e 5 'i = (fcl ^^ +fc2 ) (/W-/W), e 5 '* = (fca+ §^ +fc .) (fo«i2 - $1*22), e R * = 

(fci+fc*)(fc2+fc|) | fei+fc| p ( K n K 22-^i2K 2 i) and Kij = ^"l.+tf 3 " • The six arbitrary complex parameters ai,a 2 ,/?i,/?2,fci 
and k 2 determine the amplitude, velocity and phase of the asymptotic soliton. As we have detected already, we see also 
here that the two-soliton solution (19) of (||) differs from that of the Manakov model in phase-terms in a nontrivial 
way. Clearly, the phase change depends on the values of the real free parameters p\,p 2l T\ and t 2i and vanishes 
for the trivial choice giving back the Manakov soliton. A natural question therefore arises: Does this change in 
phase of (qi,q 2 ) for non-zero values of pi,p 2 ,Ti and t 2 , which in turn accounts for the effect of the quintic terms 
in (^|), make any qualitative change in the behaviour of soliton collision? This can be directly studied using the 
two-soliton solution (19-20) parallel to the procedure of the Manakov model given in Ref. |23| ], In order to answer 
the above question, we have plotted pictures of soliton collision corresponding to the function (|<;it| 2 , |<?2t| 2 ) instead 
of (|gi| 2 , 1 <72 1 2 ) , since due to (|<7i| 2 , 192I 2 ) = (|9im| 2 , |<7 2 m| 2 ) one has to look into the derivatives of the fields, where 
the effect of phase terms are reflected. For comparing with the pure Manakov model let us consider first the case 
pi = p 2 = n = t 2 = 0. Fig. 2, shows the asymptotic forms of (|<Zit| 2 , |<72t| 2 ) for this case at z = +7 for the parameter 
values ki = 1.5 + i0.5, k 2 = 2.0 — i0.7, o>i = Pi = /3 2 = 1, an d a 2 = ^ 39 1"q 80 ^ . At z = —7, we have two well seperated 
asymptotic profiles as shown in Fig. 2a. During the propagation, these two solitary profiles interact with each other 
as shown in Fig. 3 and change form after interaction. For example, at z = +7 they have the profiles as shown 
in Fig. 2b. The change in shape disappears if we apply the elastic collision (shape-preserving collision) condition, 
namely a\ : a 2 — j3\ : (5 2 by following the work of ]23[] . In Figs. 2-3, as mentioned, there is a splitting in each of the 
asymptotic profiles which appear before and after interaction. This happens due to the following reason. In the case 
corresponding to the Manakov model one obtains 

VjMt ~ A 7 pk nR e" M sech{j lnR + (j> nT ) { [~k nR tanh(r ]nR + </> n=F )] + ik nI ) \ z ^ ±00 ,j, n = 1,2 (21) 

where rj nR = k nR (t — 2k n iz), r\ ri i = k n it + {kf lR — k^j)z, the subscript j denotes mode while superscript n=F is used to 
define the two different interacting solitary waves appearing at z ~ =Foo and A™^ and </> n=F determine unit polarization 



vector and phase of the modes as defined in Ref. |23|. From (21) one can note that for suitable choice of the parameters 
k n R and fc„/, the solitary waves get peaked around two values as shown in Figs. 2-3. 

Now to investigate the effect of non-zero values of pi , p 2 , t\ and t 2 or in other words to see the nontrivial contributions 
due to cubic-quintic generalization (||), we evaluate qjt — [qjMt + iQjM^jt] exp(i8j), j = 1,2. and plot the 
asymptotic behavior of the two-soliton solution (19) to the quintic generalization of the Manakov model Eq.(^) in 
Figs. 4a, b. The corresponding interaction profile of the solitons during their propagation is shown in Fig. 5. We 
observe firstly that, as in the case of Manakov model, here also generically the fascinating shape changing inelastic 
collision persists. However, in this case one can overcome the splitting effect of Figs. 2-3 corresponding to the Manakov 
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model. For example, if we set p\ = pi = t\ = T2 = 1, in (19) then the splitting of solitons disappear, as evident from 
Figs. 4-5. The reason for this is that now we have 

l^i | 2 = \qjMt\ 2 + \qjM\ 2 \0]t\ 2 = \qjMt\ 2 + |<2jA/| 2 (<5l|<7lj\/| 2 +S 2 \q2M\ 2 ), 

where Si and 82 are equal to p\ and p 2 for j = 1, while to t\ and T2 for j = 2. Since here the second term dominates 
over the first in the region of splitting, the splitting effect naturally gets supressed. Comparing Figs. 3 and 5 it 
is also important to note that the intensity of solitons ( | <7x* | 2 , |<72t| 2 ) at the intersection region for solution (19) of 
our generalized model is much higher than that for (20) corresponding to the Manakov model. The above processes 
vividly demonstrate the nontrivial effect of the additional terms involving parameters pi, p 2 , t± and t 2 appearing in 
Eq.(§. 




V. RELATION WITH THE MANAKOV MODEL 

As we have mentioned above there exists an interesting interrelation between the quintic generalization (|p) and the 
Manakov model (0), which in fact we have used already in deriving the soliton solutions of (g). We establish now this 
relationship by showing that the Lax operators of these two models are related through a local gauge transformation 
p9| pl| , while the fields are connected by a nonlinear transformation in dependent variables. 

It is known |52| that under a gauge transformation of the Jost function $ = 3$ with the gauge field g <E U(3), the 
Lax operators transform as 

L' = g~ l Lg - g- l g u M' = g^Mg - g' x g z . (22) 

Choosing now the specific form 

/ 1 \ 

(23) 

with its elements 9\, 02 being the same functions of z and t as in Eq. (11c) and performing the transformation ([22]), 
one can conveniently remove the diagonal terms involving #it,#2t and 9i z ,9 2z in the Lax pair (11a, b). It can be 
observed further that the resultant gauge transformed Lax operators reduce exactly to those of the Manakov model 
p2| , [35f if we introduce transformed fields 

q a M = q a cxp(-i6> Q ), a = 1, 2 (24) 

along with their conjugates. At the same time transformation (|24|) reduces the equations (||) to those of the Manakov 
model (§). 

The above points establish the relationship between these models and justifies the form of soliton solution presented 
in the earlier section for the model (^|). Moreover, this procedure also provides an alternative proof of the integrability 
of our model. It is important to note that under such gauge transformation the Poisson bracket structure of the fields 
also gets changed. To find such changes in the canonical structure we may use transformation (|24|) to express our field 
through the Manakov fields and assuming standard canonical relation (|2^) for the Manakov model, we can derive the 
anyon-like relations for the fields of (|9|): 

{qi(x), ql(y)} = S(x -y) + ipxe(x - y)qi(x)ql(y), 

{qi(x), qi(y)} = ipie(y - x)q 1 (x)q 1 (y), 

{qi(x), 92(2/)} = -i(p 2 0(x -y)~ ndiy - x))q 1 (x)q 2 (y), 

{qi(x),q 2 (y)} = i(p20(x -y)- n6(y - x))q 1 (x)q*( y ), (25) 

etc., where e(x) = 6(x) — 0(—x) is the sign function defined through the step function: 9(x) — 1 for x > 0, 9(x) = 0, 
for x < 0. Note that at x = y the fields exhibit canonical property, while at x ^ y their behavior is nonultralocal 
and mimics anyon-like properties [34| in the classical limit. It may be remarked here that the generalized Manakov 
equation (9) can be derived directly from the Hamiltonian (17) by careful application of the PB structure (25) and 
the relation d x 9(x — y) — S(x — y). 
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VI. VECTOR GENERALIZATION OF CHIRAL SOLITONIC MODEL 



We have seen that for obtaining (||) from the Hamiltonian (|l7| ) we have to use noncanonical brackets (^j|) . On the 
other hand if nevertheless one considers them to be canonical, i.e. 

{q i (x),q*(y)}=S(x-y)S ij , { qi {x), qj (y)} = 0, (26) 



from the same Hamiltonian (17) we can derive completely different coupled equations with fifth-degree nonlinearity. 



If for simplicity we assume pi = pi = t\ = ti = po, we can derive these equations easily from C3 as 

mz + Out + 2(kl| 2 + |g 2 | 2 )gi - Sp 2 (\ qi \ 2 + \q 2 \ 2 ) 2 qi - 2ip ((\qi\ 2 + \q2\ 2 )qit + (qUlt + = (27) 

and similarly for q 2 by interchanging the indices 1 <-> 2 in (p7j). We notice that this system of coupled equations 
again with cubic-quintic nonlinearity is a new system which is different from Eq. (^) presented earlier. To anlyse 
these equations more closely we perform again a nonlinear variable change as q a — -> Q a — q a e~ lp ° e with t — N = 
\Qi\ 2 + \Q2\ 2 ■ After some lenghty but simple manipulations one can reduce the system ( p7| ) further to a more compact 
form with only cubic nonlinearity: 

iQaz + Qatt + 2(N - P0 j)Qa = 0, (28) 

where we have denoted j — j\ + j2,ja = i{Q* a Qat — QaQ* a t)- We immediately recognize that this is nothing but the 
vector generalization of the Aglietti et al equatio n [|33| , however with the addition of a cubic nonlinearity coming from 
the Manakov model. Nevertheless the system (|28|) shows remarkable property close to the chiral-soliton feature of 
|j33f . In particular assuming the 1-soliton form as Q a = A a s{t — vz)e l ^ t+u)z ^ one may conclude that here the quantity 
k = 1 + vpo acts as the effective coupling constant of the nonlinear term, which regulates the intensity of the soliton. 
Therefore for the soliton velocity v > — i only (with p > 0) Manakov model like bright solitary wave solution can 
exist. With decreasing velocity the effective coupling constant n also decreases, which interestingly influences the 
intensity of the soliton to increase and reflects a possible nonintegrable property of the model. Finally for the soliton 
velocity v = — the nonlinear term to sustain the soliton disappears and hence no such soliton can appear anymore. 

However for the negative velocity below this value, v < — the sign of k flips and kink-like exact dark soliton can 
appear. For po < the whole picture reverts. This amazing solitonic feature evidently is a generalization of the chiral 
soliton property of |33f| due to the presence of Manakov term as well as the multi-component nature and may have 
important applications in nonlinear optics. The suspected nonintegrable nature of this system and consequently the 
original chiral-solitonic system [p3| can be convincingly proved by showing that the conserved quantities of the model 



are not in involution (in particular using the canonical bracket (26) it can be shown that {02,03} ^ 0). Therefore 
though this system possesses Lax pair and infinite conserved quantities, their noninvolutiveness spoils the integrability. 
The involution of the conserved quantities however is restored if we use the noncanonical bracket ( |25| ) and this ensures 
the exact integrability of (^|). 



VII. CONCLUSION 

We have constructed the Lax pair of the proposed integrable CNLS equation (^) with cubic-quintic nonlinearity 
governing the soliton propagation in non-Kerr media, and using it generated the infinite set of its conserved quantities 
in explicit form. We also presented the exact one and two-soliton solutions of the model using those of the well 
known Manakov model. It has been demonstrated through the explicit two-soliton solution of the proposed model 
that the intensity of the t-derivative of soliton in the interaction region is much higher than that of the Manakov 
model. Moreover, the localized part of the time derivative of the Manakov soliton gets splitted and peaks around two 
values as shown in Figs 2 and 3. However, such splitting can be suppressed in the generalized cubic-quintic equation 
having nonzero p's and r's as has been demonstrated in Figs. 4 and 5. These figures also confirm that the shape 
changing inelastic soliton collisions, as in the Manakov case, persist in our model. We believe that our results will be 
found equally useful in more general situations like Eq. (5) by taking our model as the unperturbed part and treating 
the remaining terms as perturbations. 

We have also established the relationship between the proposed model and the Manakov model at the Lax pair 
level as well as at the field solution level, which shows an intriguing change in the canonical structure, namely the 
bosonic relations of Manakov model transforms into the anyonic relations of the present system. 
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Another remarkable fact is that assuming the standard canonical structure for our fields we are able to derive from 
the same Hamiltonian yet another coupled system with cubic-quintic nonlinearity. This novel model, which turns 
out to be nonintegrable, represents a vector generalization of Aglietti et al model, famous for exhibiting chiral-soliton 
solutions. Such chiral-soliton property also prevails in the present vector case showing fascinating properties of the 
solitons, like changing intensity with soliton velocity, vanishing of bright solitons and appearance of dark solitons 
below certain velocity etc. Such properties may have important applications in nonlinear optical processes. 
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FIG. 2. Asymptotic forms of the intensity profiles |giMt| 2 and |g2Mt| 2 of the two-soliton solution (20) of the Manakov model 



( pi = p-2 = T\ — T2 = 0) with the parameter values k\ = 1.5 + i0.5, k% = 2.0 — i0.7, cti = /3i = fh = 1, o<2 
z — —7 (before interaction) , (b) at z = +7 (after interaction). Note the splitting in the asymptotic soliton profiles. 
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FIG. 3. Intensity profiles |gu\/t| 2 and \q-2Mt\ 2 of the two-soliton solution of the Manakov model with the parameteric values 
as in Fig. 
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FIG. 4. Asymptotic forms of the intensity profiles \qit\ 2 and \q2tf of the two-soliton solution (19) of the generalized model 
(^), (a) at z — —7, (b) 2 = +7 , for nonzero values of the parameters pi = p2 = r\ = r% = 1 and with remaining parameters as 
in Fig. 3. Note the suppression of the soliton splitting, which appeared in the asymptotic profile in Figs. ^ and |^. 
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FIG. 5. Intensity profile \qit\ 2 and \q2t\ 2 of the two-soliton solution (19) of the generalied Manakov model with the parameter 
values as in Fig. 4. Note the persistence of inelastic soliton collision as in the Manakov model and a higher intensity of modes 
during soliton interaction compared to the Manakov model. 
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